Mixed scenario of the charged 
helium surface reconstruction 



V.Shikin, E.Klinovaya 
Institute of Solid State Physics of RAS, Chernogolovka, Moscow distr., 142432, Russia 

Abstract 

Discussed in the paper is the mixed scenario of charged hquid surface recon- 
struction when the surface 2D charge density is close to the saturation. The basic 
building block of arising honeycomb structure is shown to be a modified multielec- 
tron dimple. 
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One of the instabilities studied in classical hydrodynamics is the the Frenkel-Tonks 
(FT) instability [H |2] occurring in a threshold manner at the charged liquid surface and 
resulting in its deformation. A distinctive feature of the FT instability compared with 
other decay processes: (Rayleigh instability of the cylindrical jet [3], Karman trace behind 
the moving cylinder (sphere) instability of the interface between two liquid media 
moving relative to each other [5], Taylor vortex instability of viscous liquids between 
two rotating coaxial cylinders [5], etc. is the possibility of stopping the decay process 
after which a new metastable state with a finite corrugation amplitude is formed (i.e., 
a reconstruction occurs). At different times this phenomenon was considered by many 
authors for dielectric and charged interfaces [6] - [28]. However, many important details 
of this process have only benn understood quite recently. In the FT instability, only the 
initial part of the problem (calculation of the dispersion law and it stability analysis) 
has a universal nature. Further analysis is requires different treatments for the dielectric 
and charged boundaries. In contrast to the dielectric interfaces, the latter case allows 
formation of the solitons of a special kind, the multielctron dimples, [211 [22] . A less 
obvious, but actually more important factor resulting in different behaviour of nuetral 
and charged surfaces is the violation of homogeneity of boundary conditions (the surface 
charge distribution along the interface can become discontinuous). Earlier Mel'nikov and 
Meshkov showed ^23j, i25j that the boundary condition inhomogeneity inevitably occurs in 
the problem of the charged helium surface reconstruction. Therefore, its analysis should 
account for arising inhomogeneities. 

The first of all inhomogeneous reconstruction scenarios we are aware of substantially 
employing the elements of equipotential reconstruction theory [TT], |T2l [131 [13 l2Qj, and 
hence the important for this formalism low electron population z/ <C 1 of the liquid 
surface was proposed by authors of Ref. [231 125]. Omitting the details, we only note 
that "bald" areas arise in their treatment as a triangular lattice (whose primitive vectors 
are equal to the capillary lengths a) of small spots having the shape of circles of radius 
R* <^ a, a = K~^, growing from zero as the supercriticality rises. Unfortunately, the 
theory given in Refs. [221 [25] does not contain individual solitons (which are typical of 
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the charged hehum surface reconstruction observed for z/ <C 1) among its solutions. In 
addition, the maximum surface curvature in the central parts of nuetral spots has the sign 
opposite to that observed experimentally [27]. These circumstances hamper evaluation 
of the feasibility of the model proposed in Ref . [231 125] . Later, it was found that in the 
range of z/ ^ 1 the details of the reconstruction are much more naturally described in 
terms of individual multielectron dimples [28]. Here and below the boundary population 
with mobile charges v is defined as the ratio 

z/ = n,/r^^^ 27r (enr")^ = (1) 

where p and a are the density and surface tension, g is the acceleration due to gravity, e 
is the electron charge, k = ^Jpg/a is the liquid capillary constant. 

In the opposite, most suitable for observations limit u < 1 where the equipotential 
language [151 120] and its modification [231 125] do not work at all, it was suggested to 
model the corrugation with inhomogeneous electrostatic potential distribution by a pe- 
riodic system of multielectron dimples [26]. The analysis of Ref. [26] is based on tight 
binding approximation. The structure is assumed to consist of roughly estimated free 
multielectron dimples. By further assuming these quasiparticles to be point-like, the au- 
thors model the periodic lattice of charged dimples by employing the results derived for 
classical Coulomb crystal consisting of point charges [29]. The dimple crystal spacing is 
approximately equal to the capillary length while its advantage compared with the homo- 
geneous 2D charge distribution is treated in terms of correlations [29j. To comment on the 
results of Ref. [26j, we note that the single dimple problem accurately solved in Ref. [21] 
yields for the conditions of Ref. p6] the charged spot radius considerably exceeding the 
lattice spacing given in Ref. [26] (according to Ref. [21], here R/a ^ 1.5). The critical 
conditions for its development is higher than the dynamic stability threshold. Finally, 
we believe the correlation effects [29] are not directly related to the Coulomb part of 
the reconstruction problem. Thus, the dimple picture of a periodic reconstruction, being 
almost the only realistic in the range of u < 1, contains in the formulation given in Ref. 
[26] a number of disputable concepts whose presence should be at least understood and 
which should be modified. The details of such an analysis are given in the present report. 

1. From our point of view, most questionable are suggestions of authors of Ref. [26] 
on the Coulomb part of the problem. It is assumed that the development of a periodic 
corrugation is strongly assisted by the well-known gain in energy occurring when the 
2D charged systems goes from the gaseous (liquid) state to the crystalline one (Wigner 
crystal). Such a gain of the correlation origin indeed takes place (e.g. see Ref. p9]) as 
revealed by analysis of the difference between the average Coulomb energy of (i) likely 
charged point-like particles < V^{r) > uniformly distributed along the surface with the 
average density Ug, and (ii) those arranged in a lattice < V^{r) > with the same average 
density 

S < Vcir) >=< 1/'(r) > - < r/(r) >~ +0,6^7^. (2) 

Here the constant Cg of the order of unity depends on the lattice type. By employing a 
formal analogy between the structure of a classical Wigner solid with one electron per cell 
and the corrugation having spacing b and the charge go per cell, it is possible to assign to 
the Coulomb part of the reconstruction problem (as it was done in Ref. [26]) the energy 
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Figure 1: Energy gain SW ([7]) as a function of the electron spot radius. The negative 
minimum corresponds to the optimal value of R/a. 



gain of the order of 

6 < W,ir) >~ +Csq^o/b > (3) 

However, the hypothesis Q can actually be checked and, at least for the employed 
model of a corrugation involving formation of a set of one-dimensional conducting threads 
in the practically interest case u = 1 is not confirmed. Detailed calculations (which can 
be found in the Appendix) lead to the energy whose sign is opposite to that in equation 
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S < W,{r) >= ^ (l/C - 1/C) < (4) 

~ And' ~ ATTd + Aa\n{a/TTR) 

where S is the capacitor area, R is the single thread radius, 2a is the corrugation period. 
In addition to the sign which is opposite to that of Eq. [3]), the formula (jl]) has one more 
qualitatively important property: here the energy gain depends explicitly on the length R 
measuring the charges localization length in a separate groove. This length is determined 
self-consistently by requiring the minimum of the corrugated system total energy. In the 
approximation used in Ref. ([3]), the Coulomb part of the problem which does not contain 
R slips out of the self-consistency condition which should not be allowed. In the rest of 
the paper we follow the definition of 6 < Wc{r) > given by Eq. (jlj). 

2. Here are some remarks before the study of the corrugation problem for u < 1. First 
of all, the total charge Q at the liquid-vapor interface with area S* ~ is constant: 

Q = const = enr'^L^ (6) 

Perhaps a trivial point, but the problem formulation should for definiteness contain 
some comments on the initial stage of the decay. It develops within the dynamic equipo- 
tential scenario [TB], until under the conditions of Eq. the growth of the corrugation 
amplitude results in the break of continuity in the charge distribution. In the limit of 
small u the inevitable appearance of continuity breaks was noted by authors of Refs. 
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Figure 2: Fig. 2. One-dimensional dimple profiles for a single free dimple (dashes) and a 
periodic set of dimples (solid line). 



[231 [23] . However, their reasoning is essentially universal in nature and is applicable for 
arbitrary values of v. If one assumes that the spatial period of developing corrugation is 
not changed with time (in the general theory of spinodal kinetics relevant to the problem 
considered here this assumption is not true [30]), the net charge splits into a system of N 
individual clusters of charge Qo per unit length 

iV = L/2a, Qo ^ 2aenr'' (7) 

each of which is the nucleus of a separate one- dimensional (for simplicity) dimple residing 
at its center. It should be noted that the unstable state could decay in a purely fluc- 
tuational (binodal) way just as in the case of z/ <^ 1 (see Ref. [28]). However, in the 
general case of arbitrary v the dimple charge remains undefined in this scenario. The 
decay having a spinodal initial and binodal final stages (we call it "intermediate") does 
not contain this uncertainty. 

In the problem of finding the corrugated surface state in the critical electric field E"^"-^ 
with the charge Qo ^ per cell, one should verify the ability of the deformed helium 
surface to hold the critical charge (perhaps tking into account the non-linear phenomena) 
and check the validity of inequality 

iy(E!;^"^ Q) > ^y(E!;^'^^ Qo, &), ^r"" = 47renr^ (8) 

where n™"^ is taken from Eq. ([1]), WiE"^"-^ ^Q)^ W [E'^"-^ ^Q^^h) are the total energies 
of the uniform and corrugated helium surface states, h is the lattice spacing which is 
generally different from the capillary length a. There are no reasons to require coincidence 
of energies in Eq. ([H]) since we do not consider the actual binodal transition from a flat 
to corrugated state. Nevertheless, it is desirable to have the energy W^E"^"-^, Qq, h) to be 
less than W{E^"'^,Q), since otherwise the geberal picture of the system evolution could 
not be reasonably interpreted. 

Bearing in mind the above arguments and results derived in the Appendix, the problem 
of one-dimensional corrugation is reduced to calculating the energy difference W{E'^°'^, Q) - 
W{E^°-^, Qo,b) ([8]). Note that the one-dimensional scenario is not only a convenient model 
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but also an experimentally observed state of the corrugated surface [27] . In addition, this 
is the least stable type of corrugation. Laplace pressure stabilizing the picture is twice as 
low in the one-dimensional case as in two dimensions. The existence of a solution yields 
an estimate from below for the corrugation stability. 

The initial energy functional for a one-dimensional corrugation (per unit thread length) 



-l-u 

5Ty (E!^'^^ i?, a) = J dx 



f[(V0' + «'e']+ei^;r'5n(a;)e(x) 



6 < W,{r) >= 2e%y In 



+ 5<W,ir)>, (9) 



^(x) = ({x)— < ( >, 6n{x) = n{x) — Ug 

represents the system energy counted from the energy of uniformly charged fiat liquid 
surface and written for one lattice cell. Its Coulomb part is taken in the form of Eq. (4) 
allowing variational solution for the problem of finding the length R if the charge density 
n{x) in a periodic system of dimples is split into individual stripes each having some 
density distribution no(x) (for example, a Gaussian one) 

n{x) = y2^o{x -2la), no(x) = no exp (-xV^?^), V^^rio = — „, ' , ^ (10) 
^ erj[a/n)K 

(erf(x) is the Icnown special function) normalized according to Eq. (I7l). Average values 
of the variables C(^); i^lx) are subtracted from their initial definitions according to the 
general structure of the functional ([9]). 

The functional (Q has the standard for the theory of multielectron dimples structure. 
The deformation-related part written in square brackets makes the reduction of the pa- 
rameter R/a as well as the growth of liquid surface deformation due to local electron 
pressure energetically favorable. The Coulomb component 6 < Wc{r) > hinders this pro- 
cess. The interplay between these factors results in the appearance of a negative minimum 
in the 6W{E'^"'^, R/a) (Q dependence which determines the equilibrium value of the pa- 
rameter R/a. To get the numerical values, one should solve the liquid surface mechanical 
equilibrium equation for ^(x) assuming n{x) (fTOj) to be a known function, then calculate 
the integrals with these distribution, and finally obtain the energy dependence on 
the variational parameter R/a. For the external parameters corresponding to u = 1 this 
picture is presented in Fig. [H 

The functional minimum is reached at R/a = 0.1 and is negative. In other words, 
the corrugated state is more favorable than the flat one. The corrugation is stabilized in 
the range of R/a <C 1, although the obtained numbers should only be considered as a 
rough estimate because the values of (x) < 1 are not parametrically small, as follows 
from the data on the one-dimensional dimple profile presented in Fig. [2]). Also shown for 
comparison in the same figure is the individual free dimple profile. Obviously, the linear 
theory corresponding to the bilinear energy functional (Q requires quantitative corrections 
(it follows from the figure that the inequality V^(x) < 1 used to derive Eq. ([9]) is only 
poorly met). However, qualitatively the mixed reconstruction scenario is self-consistent 
(starting from the assumption on the presence of continuity breaks in the equilibrium 
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distribution n{x) we finally obtain the picture confirming the initial assumptions of the 
theory) . 

To sum up, we conclude that the mixed scenario of the charged liquid surface recon- 
struction under the conditions of u = 1 involving the formation of a periodic system of 
charged dimples has a self-consistent nature. The period of arising structure is determined 
by the dynamical instability of charged liquid surface. Development of this instability di- 
vides the surface charges into a periodic system of clusters localized in every cell of a 
one-dimensional or a honey-comb corrugation. When approaching the metastable state, 
the charged cell cores becomes rather small {R/a -C 1) while the corrugation amplitudes 
is of the order of the capillary length (stiff reconstruction regime according to the clas- 
sification of Ref. [H]). The mixed scenario reveals no hysteresis phenomena typical of 
nonlinear equipotential theories (which is not surprising because of the employed linear 
approximation (9) for energy). 

Formally, the Coulomb energy should be treated in a consistent way. The field- 
theoretical approach allows one to avoid paradoxes at both quantitative and qualitative 
levels. In particular, in the considered problem there are no serious reasons for manifesta- 
tion of correlation phenomena governing transition to crystalline ordering of the Wigner 
type. 

This work was partly supported by the RFBR grant 09-02-00894. 
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Appendix 

The statements from Ref. [28] concerning the properties of the difference 

S < V,{r) >=< V^ir) >-< V;(r) > (1) 
are based on the analysis of the expression 



6 < Vc{r) > = lim ^ 
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\x — x{l)\ \x\\ 
x{l) = hai + ^202, (3) 

where ai 02 are the primitive translation vectors in the considered lattice. 

The sums in Eq. ([2]) diverge already because the classical proper energy of every 
electron in the definition ([2]) is infinite. However, the Fourier analysis and the important 
condition of the system neutrality as a whole (which is insured by introducing the screening 
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electrode) allow to reduce the difference of the sums to a finite expression of the form 
(e.g., see Ref. [22]) 

S < K(r) >= Ncse^y^j2, (4) 

where is the total numcer of elctrons in the 2D system, Oc — is the unit cell area, 
Cs is a constant of the order of unity depending on the lattice type. 

An elegant and rather general result (jlj) gives an idea on the nature of Coulomb 
crystallization of electrons above helium yielding the crystal ground state energy in this 
problem compared to the gaseous (liquid) phase. 

Let us now turn to the known electrostatic problem [31] on the field and energy of a 
set of charged threads above a conducting screen, i.e. the problem of a fiat capacitor with 
a solid or mesh-like electrodes. Here the electrostatic potential has the form 

, , ^ , f sin In (z — id)/ a\\ . 
^{z) = -2aalni . \ ) ^ .J' , z = x + ly 5 
\_sm [7r(z + id)la\ j 

where a is the lattice spacing, a is the mesh charge per unit area, d is the distance 
between the mesh and the screen. Just as in ([2]), the divergencies at infinity cancel each 
other. The filed is concentrated within the capacitor, the mesh filed becoming uniform 
at distances of the order of its period 2a -C Id where Id is the distance between the 
plates. The total electrostatic energy (the counterpart of Eq. ([2|)) calculated in the field- 
theoretical approach (squared electric field integrated over the capacitor volume) proves 
to be divergent due to the potential i{>{z) ([5]) singularity at the threads. This property 
requiring cutoff at the thread radius R does not appear in the final result (jlj). 
Capacities C , C per unit system area after using Eq. (jS]) can be calculated to be 

c = ^, c ' 



47r6' 47rb + 2aln{a/27rR)' 

Here i? <^ a/27r is the radius of a single thread. 
The counterpart of expression ([2]) is obvious: 

6<W,{r)>=^il/C-l/C) (6) 

This difference has the sign opposite to that of Eq. (j4]) and depends on particular 
characteristics of a single thread. Omitting any comments on possible reasons of the 
qualitative difference between Eqs. (jlj) and ([6]) we shall employ in the main body of the 
paper the Coulomb part defined by Eq. ([6]) since it contains information on the internal 
Coulomb energy of every charged thread which is absent from Eq. (jl]). 
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